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Random quantum spin systems offer a useful laboratory for studying the fascinating 
interplay between strong interactions and disorder. Though not as complex or intractable 
as metal-insulator transition systems, they are still rich enough to display a host of unusual 
physical phenomena. Moreover, they can be realized in a number of experimental systems, 
many of which have been studied intensively in recent years []!]-§[. 

It is useful to distinguish two different types of possible ground states of a random 
quantum magnet: (a) a state with magnetic long-range-order ((Si) ^ where <Sj is the 
spin operator on site i) which can be a spin-glass, ferromagnet or an antiferromagnet; (b) a 
quantum disordered (or 'spin-fluid') state in which (Si) = due to the presence of strong 
quantum fluctuations. Many properties of the magnetically ordered phase can be described 
by a semi classical analysis. In contrast, the spin-fluid phase and its zero-temperature phase 
transition to the magnetically ordered phase are intrinsincally quantum mechanical, and 
their properties are only very poorly understood. This paper shall mainly focus on the 
properties of the spin-fluid phase. 

In early studies of random-exchange spin- 1/2 Heisenberg spin chains by a numerical 
renormalization group method, Ma et.al. and others [[J noted that the low temperature 
spin susceptibility x(T) behaved approximately like T~ a with a < 1. This behavior, and 
their analysis, suggested that the quantum disordered phase generically possessed gapless 
excitations: the low energy excitations arose from a significant probability of finding a pair 
of spins which were essentially decoupled from the rest of the system, and with only a weak, 
mutual, effective exchange interaction. Subsequently, the results of a numerical analysis by 
Bhatt and Lee [|7| of a three-dimensional random-exchange spin- 1/2 Heisenberg antiferro- 
magnet could be well fit with the same functional form with a ~ 0.66. Experiments || 
on many lightly doped semiconductors have also found similar behavior in the low temper- 
ature spin susceptibility. More recently, Doty and Fisher |J|IIJ have obtained numerous 
exact results on random quantum spin chains; in particular, Fisher JTD] proved that the 
random-exchange, spin- 1/2 Heisenberg chain has x ~ ^-/(T\og 2 (l/T)) and is gapless. 

In this paper we introduce a new solvable, random-exchange, quantum Heisenberg mag- 
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net - its solution reduces to the determination of the properties of an integro-differential 
equation, which is a difficult, though not impossible task. Our model possesses infinite- 
range exchange interactions, and is thus a solvable limit which is complementary to the 
spin chains. The spin-fluid phase of our model is generically gapless; however the physical 
mechanism of the gaplessness appears to be quite different from that of the random spin 
chains. Which of these two limits is closer to realistic three dimensional models remains an 
open question. Finally, our model is expected to display a transition to a spin glass phase. 
We have not yet succeeded in unraveling the nature of this transition and that of the replica 
symmetry breaking in the spin-glass phase - these are issues we hope to address in a future 
publication. 

The main result discussed in this paper is that the T = 0, average, local dynamic spin 
susceptibility of our model has the following form over the entire quantum disordered phase: 



l0g ( Ju\ ) +i f sgn ^) 



+ ••• (1) 



where X is a constant to be determined below, and the omitted terms are subdominant in 
the limit \u\ —>■ 0. A notable feature of this form is that it is identical to the 'marginal' 
Fermi liquid susceptibility proposed on phenomenological grounds by Varma et. al. [II] as a 



description of the electronic properties of the cuprates. It is not completely unreasonable to 
begin a study of the low-lying spin fluctuations in the cuprates by using the infinite-range 
quantum spin model described below; however, at present we have no arguments which can 
determine whether, or how, the 'marginal' spectrum will survive in more realistic models 
with charge carriers and finite-range interactions. Nevertheless, to our knowledge, ours is 
so far the only bulk model to display the 'marginal' spectrum over an entire phase, and one 
might hope that mathematical structure of the mean-field theory is of broader significance. 
We consider the ensemble of Hamiltonians 

where the sum over i, j extends over N sites, the exchange constants J^- are mutually 



uncorrelated and selected with probability P(Jij) ~ exp(— J^/(2J 2 )), the <S are the spin- 
operators of the group SU (M), and the states on each site belong to a representation labeled 
by the integer rib {rib = 2S for SU(2)\ more generally rib is the number of columns in the 
the Young tableau of the representation [EL3J). This model has been considered previously 



by Bray and Moore (TJJ for the group SU(2): they found strong evidence in favor of the 
presence of spin-glass order at T = for all values of S. Accessing the spin-fluid phase 
therefore requires considerations of groups other than SU(2): following a technique which 
has been successful in clean antiferromagnets [fT3],|l4T] , we generalize to the group SU(M) 
and study the phase diagram in the rib — M plane. (We have also studied the properties of 
random Sp(M) Jl5| magnets with results that are very similar to the simpler SU(M) case 
considered here). There are three interesting limits in the rib — M plane: 

(A) rib —>■ oo, M fixed: this is the semiclassical limit and yields ground states well within 
the magnetically ordered spin-glass phase. 

(B) M — > oo, rib fixed: this takes us deep into the spin-fluid phase. 

(C) M — > oo, rib/M = k fixed: this is in many ways the most interesting limit, because by 
varying k one can interpolate between the spin-glass and spin-fluid phases. Moreover, one 
expects a phase-transition between these two ground states at a critical value of k — k c . 

The structure of the N —>■ oo limit was discussed in Ref. [0. We express the partition 



function as a coherent-state path-integral [13], introduce n replicas, average the partition 



function, and the saddle-point reduces to the quantum mechanics of n replicas of a single 
spin; assuming the saddle-point is spin-rotation invariant (this is true in both the spin- 
fluid and spin-glass phases) we obtain the single-site coherent-state path integral Z = 
J VS exp(£) with 

C = S B + / o VT drdr'Q ab (r - r')S a (r) ■ S b (r') (3) 
and the self-consistency condition 

Q ab (r-r') = ^(S a (r).S b (r')) Zo (4) 



Here a,b = 1 . . . n are replica indices, r, r' are Matsubara times, and Sb is the single-spin 
kinematic Berry phase term The Edwards- Anderson order parameter [16[] for the spin- 



glass phase is g^M = Q aa (r —>■ oo). Moreover, Q ab , a 7^ 6, is r-independent and non-zero 
only in the spin-glass phase fT/TJIBfl . 

An exact evaluation of Z Q is clearly not possible. We therefore consider the large M limit, 
discussing first the limit (C) above. This is achieved by the Schwinger boson realization of 
S 

S a ; = b a X ; ]T6^ = n 6 (5) 

n 

where h is a boson annhilation operator, //, v = 1 . . . M. In the large M limit, 
Eqns (|3]|D reduce to the following equations for the boson Green's function G a B {r) = 



(1/M) Eu(T(b atl {TW(0))) and its Fourier transform G 



ab 1 



G B {iu n ) = (-iuj n + A - E B (iu n )) 1 (6) 



£#(r) = J 2 G a B b (r)G a B b (r)G b B a (-r) (7) 

while Q ab (r) = G a B {r)G b B {—r). Here A is a chemical potential set by the constraint 
G aa (r = 0~) = k. These two equations can be combined into a a single integro-differential 
equation for G a B {r). We also require that solutions satisfy conditions imposed by the 
spectral representation of a boson Green's function: G B a (z) is analytic for Im(z) > 0, 
uj lm(G < B (u + i0 + )) > and G a B (z) ~ — 1/z for large \z\. The replica-diagonal compo- 
nents of Eqns (|||7|) also bear a formal resemblance to a perturbative solution of the infinite- 



dimensional Hubbard model |19[ ; however there are some significant differences which turn 
out to have dramatic consequences in the nature of the solution. 

We will focus here only on the spin-fluid phase, whence all correlations are replica- 
diagonal, and replica indices will be dropped. An immediate consequence of (||J^) is that 
the zero-temperature boson spectrum must be gapless! For suppose that the spectral weight 
Iul(Gb(u; + i0 + )) = for |u;| < A; then (^), expressed in real frequencies, implies that 



lm(E B ((j + iO + )) = for \u\ < 3 A - this agrees with the real- frequency version of (||) only 
if A = 0. 

Let us focus on the low- frequency behavior of Gb'- assume that Gb{oj) ~ u/\ Then from 
(0) we get Im(£#(u;)) ~ u 2+ ^. This can be consistent with (^) only if A = £_b(cg> = 0) and 
/i = —1/2. As Gb is analytic in the upper-half frequency plane, we write 

G B {z) = —=- + ... lm(z) > (8) 

V z 

where A > 0. The positivity conditions on the spectral weight require < 9 < n/2. Inserting 
this into (0) we find for Im(z) > that 

^ / x _. > «. J 2 A 3 sin(2#) m r - 

7T 

Finally, this is consistent with (Q) if A = S^(0) and 

/ \ 1/4 

The parameter 9 remains undetermined. This is fortunate, as we need a single degree of 
freedom to satisfy the boson-number constraint Gb{t = 0~) = n. We will treat 9 as the 
independent parameter, with k{9) a function to be determined. We expect k — > 0, as 
9 — ► 0; increasing 9 therefore corresponds to increasing 'spin'. We can also determine the 
low-frequency behavior of the spin-susceptibility x(r) = Q aa {r)\ we find that it has the form 
(H) with the constant X given by 

= («g 

2J V ; 

We expect the low-frequency susceptibility to increase monotonically with increasing 'spin' 
k, and therefore increasing 9. However, X has a maximum at 9 = tt/4. This leads us 
to conjecture that the transition to the spin-glass phase occurs at 9 = tt/4 and only the 
range of values < 9 < 7r/4 correspond to the spin-fluid phase. A second possibility, which 
cannot be ruled out, is that there is a first-order transition to a spin-glass phase at a value 
of 9 < tt/4. 
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We have performed a detailed numerical study of Eqns (|6],|7|) to determine the complete 
frequency dependence of Green's function. We chose a value of 9, and a trial form for 
Iid.(Gb(oo + i0 + )) whose low-frequency limit satisfies Eqn (|8|). The real-frequency version 
of (0) expresses Im(Ss(ci; + i0 + )) as a double convolution of Im(Gs(cj + i0 + )); these con- 
volutions were performed by direct numerical integration. The real part Re(E(a; + i0 + )) 
was determined by a Kramers-Kronig transform, and A was set at A = = 0). Finally 

Iiii(Gb(w + i0 + )) was determined from and the whole procedure was iterated, until the 
solution converged. The singularities in Gb and at low frequencies were accounted for 
by performing the numerical integration in a variable x ~ a/cJ at the integration end-points 
- this absorbed the leading singularity. Subleading singularities were treated by using a 
dual mesh-size in the integration - a very fine mesh (x spacing = O.OOOSv^J) was used at 
the end-points and a coarse mesh elsewhere. Upto 1700 points were used in the numerical 
integration. There was little difficulty in converging to a solution for values of 9 less than 
approximately 7r/6; we are reasonably certain that there are physically sensible solutions of 
([|[7p for this range values of 9. One such solution, at 9 = tt/12 is shown in Fig [I] which 
was found to have k = 0.051. The numerical iteration became increasingly unstable with 
increasing 9 and did not converge to any smooth solution for large 9. Our numerical ex- 
perience is consistent with the conjecture that there are no physically sensible solutions for 
9 > 7r/4 - this is the range of values of 9 where we expect a spin glass phase. 

A complementary picture of the spin-fluid phase can be obtained by studying the large 
M limit (B). This takes M — » oo at fixed spin - one is then in a region of the phase 
diagram well away from the transition to the spin-glass phase. For technical reasons it is 
also necessary to introduce of order M rows in the Young tableau of the spin representation; 
this is discussed in some detail in Ref. ||13| . We will focus on the particle-hole symmetric 
representations which have n& columns and M/2 rows as realized by the following operator 
decomposition 

s™ = Y.f«U aau ; E/S/* = » (12) 
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where / is a fermion annhilation operator. The fermions carry replica, spin, and 'color' 
indices a, (3 = 1 . . . n^. The subsequent analysis parallels closely that for the bosons. The 
fermion Green's function is replica and color diagonal and its only non-zero component is 
G f (t) = (1/M) En(T{f aafl {T)fZl (0))). The only changes are that Eqn (fl) is replaced by 

E F (r) = -J 2 n b G 2 F (r)G F (-r) (13) 

and the positivity constraints on the fermion spectral weight is Im(Gj?(w + i0 + )) > 0. The 
presence of particle-hole symmetry requires that 1to.(Gf(u + i0 + )) is an even function of u - 
this simplifies the analysis considerably. The low-frequency limit of Gp can be determined 
completely: 

/ n \ 1/4 (_]. _|_ ^ 



'Z 

The dynamical susceptibility is found to have the same low-frequency dependence as in ([[]), 
with the constant X now given by 



X = ^ (14) 

As expected, X is a monotonically increasing function of n^. A complete solution was 
obtained numerically and the results are shown in Fig [2| 

The key unresolved issue in this work is of course the range of validity of the dynamic 
susceptibility in Eqn ([!]) - this is important in determining the significance, if any, of our 
results for dynamic neutron scattering experiments on random antiferromagnets |1|-|3[]: (a) 
What are the consequences of 1/M fluctuations in the infinite-range model Z (Eqns (0,|)) ? 
This question has been answered for a simpler infinite-range quantum spin-glass [H| where 
it was found that 1/M corrections did not modify the low-frequency behavior of the spectral 
weight. The structure of the fluctuations about the present mean-field theory is much more 
involved, but it is reasonable to expect that a similar phenomenon will occur here, (b) 
Is there an upper-critical dimension above which the properties of Zq describe the spin- 
fluid phase or its phase transition to spin-glass order in antiferromagnets with finite-range 
interactions ? (c) How are these results modified in ensembles with a nonzero average Jjj ? 
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FIGURES 

FIG. 1. Spectral weights of Gb and \ f° r the bosonic representation (|5|) at 9 = n/12. The 
'spin' is rib is order M: rib = kM, and for this value of 9, we found k = 0.051. The sum-rule for \ 
is / °° duilm(x(uj)) = vtk(1 + k) 



FIG. 2. Spectral weights of Gf and x f° r the fermionic representation (12). Now the spin rib 
is of order unity, and upto rescaling, the solution has the same form for all rib- The sum-rule for x 
is / °° dujlm(x(uj)) = nrib/4: 
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